THE DISTRIBUTION OF GEODESIC EXCURSIONS 

INTO THE NEIGHBORHOOD OF A CONE 
SINGULARITY ON A HYPERBOLIC 2-ORBIFOLD 

ANDREW HAAS 

Abstract. A generic geodesic on a finite area, hyperbolic 2-orbifold 
exhibits an infinite sequence of penetrations into a neighborhood 
of a cone singularity of order k > 3, so that the sequence of depths 
of maximal penetration has a limiting distribution. The distribu- 
tion function is the same for all such surfaces and is described by 
a fairly simple formula. 



1. INTRODUCTION 

Many of the concepts of number theory can be realized in a geometric 
setting. This is especially true in the field of classical diophantine 
approximation, where much of the theory has been reformulated as 
statements about the depth of penetration of a geodesic into the cusp, 
or non-compact, end of the Modular Surface. 

The Theorem of Bosma, Jager and Weidijk H, called the Lenstra 
Conjecture, describes the distribution of values 9 n = q n \q n x — p n \, for 
almost all real numbers x, where p n /c[n is the n th continued fraction 
convergent to x. Bosma proved a closely related result for approxi- 
mation by mediants jSJ. In [S] we proved analogues of these results 
which describe the distribution of the depth of maximal penetration 
of a generic geodesic into the cusp end of a finite area hyperbolic 2- 
orbifold. These purely geometric results were then applied to prove 
Fuchsian group versions of JH] and j^j- The distribution of the 9 n , de- 
fined with respect to a finite area Fuchsian group, is given independent 
of the group and agrees with the classical result. 

In this paper we look at how a generic geodesic on an orbifold surface 
"approximates" a cone singularity of the surface. This type of investi- 
gation was begun in [Hj. In that paper, the geometric Markoff theory 
was generalized to this setting. Here it is shown that the sequence of 
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depths of maximal penetration of a generic geodesic into a neighbor- 
hood of a cone singularity has a limiting density which is explicitly 
computable and the same for all finite area surfaces and most orbifold 
surfaces. As in jH], the resulting densities bear a striking similarity to 
the ones found in [S| and j^j, although they are somewhat more com- 
plex and arcane. The similarity becomes particularly clear when the 
distributions are expressed in terms of area. 

1.1. The main results. We shall begin with a heuristic treatment of 
some of the definitions and then state the main results of the paper. A 
hyperbolic 2-orbifold S is represented as the quotient of the hyperbolic 
plane H by a Fuchsian group G. We shall assume that S has finite area 
and, if S is compact, G is not a triangle group. A point p on S is a 
cone point of order k > 3 if, under the projection from H to S , p is 
covered by a point p which has non-trivial stabilizer of order k in G. 

An r-excursion e, of a geodesic ray 7 into a neighborhood of the 
cone point, is an arc of 7 associated to a point of maximal penetra- 
tion of the ray 7 into the radius r neighborhood, B r (p), of the cone 
point. The depth of the excursion e, denoted d(e), is the "distance" 
from e to p . When r is sufficiently small, the excursions are simply the 
arcs 7 n B r (p) and d(e) is the actual distance. The precise definition 
becomes more complicated with larger values of r. An r-excursion e, 
which is associated with a closed arc of 7 that loops about p, is called 
an approximating r-excursion. e is called an approximating excursion 
if it is an approximating r-excursion for some value r. These definitions 
parallel the geometric analogues of mediant and continued fraction ap- 
proximations, as they appear on the classical Modular surface. 

The r-excursions of a ray 7 naturally form a sequence {e^}, which 
goes out the end of 7 if the sequence is infinite. In a similar fashion 
the approximating excursions of 7 form a sequence {e*}. For a generic 
ray 7, both sequences are infinite for all r > 0. 

Let j^A denote the cardinality of the set A. For r>0, < z < r and 
a geodesic ray 7 with infinitely many r-excursions ej into the radius 
r neighborhood of a cone point p of order k, the distribution of the 
depths of r-excursions is defined as 

(1) dist k (r, z)(j) = lim \l<j<n, d(e j ) < z}, 

n— too fl 

if the limit exists. Similarly, if e* is an infinite sequence of approximat- 
ing excursions of 7 into the neighborhood of a cone point p of order k, 
then the distribution, dist^z)^), of the depths of approximating ex- 
cursions is defined just as the function ((H) above, with e* replacing ej in 
the definition, and < z < r k , where r k = sinh _1 (cot |). Note that the 



The distribution of geodesic excursions 



3 



function does not depend on a parameter r. In Corollary El of Section 
14.2.21 we show that every approximating excursion is an r-excursion 
for some r < r k . 

We are now in a position to state the main results of the paper. 

Theorem 1. For almost all geodesic rays 7 , for r > and < z < r 
the distribution dist k {r, converges to a function dist k {r,z). When 
r < r k , 

sinh z 
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The result for approximating excursions is given in 
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Theorem 2. For almost all geodesic rays 7 and for < z < r k the 

distribution dist* k (z){pf) converges to a function dist* k (z), which is given 
by the formula 
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1.2. Outline of the paper. There are two main pieces to the proofs 
of Theorems ^ and 121 First in Theorem El we prove the existence of the 
distributions and describe their values by certain integrals. In Theorem 
01 we complete the proof by computing the values of these integrals. 

Section |21 is concerned with background material and a precise treat- 
ment of excursions. The statement of Theorem |3] is in Section 13.21 and 
the proof occupies the remainder of Section El Theorem 0] is the focus 
of Section HI In the last section, Section we see how the results look 
if area, rather than distance, is used to define the depth of an excur- 
sion. Written in terms of area, it is easier to see how the work with 
cone points relates to the more classical case. 
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2. Excursions into a cone neighborhood 

2.1. Basic notions and definitions. The quotient of the Poincare 
upper-half plane HI by a Fuchsian group G is called a hyperbolic 2- 
orbifold. Let n : H — ► M/G = S denote the projection from the 
upper-half plane to the orbifold S. A point p e S is called a cone pomt 
o/ order k if there is a primage p of p in H so that the stabilizer of p in 
G is generated by a transformation of order k. 

Throughout this paper we shall assume that S is a finite area hy- 
perbolic 2-orbifold with a cone point p of order k > 3 and that, if S is 
compact then G is not a triangle group. One important consequence 
of this last hypothesis is that, if S is compact, then there is a simple 
closed geodesic on S that does not pass through p . This follows from 
the fact that such a surface will contain a homotopically non-trivial, 
simple loop in the complement of the cone points on S , which does 
not go around a cone point We should mention that in the special case 
where G is a group with signature (0; 2, 2, 2, k), the simple closed geo- 
desic degenerates to a geodesic ray going back and forth between order 
two cone points. 

For r > 0, let B r (q) denote the open hyperbolic disc in EI of radius 
r centered at q. When r = set B (q) = {q}. For p covering the cone 
point p , B r (p) projects to what we shall call the cone neighborhood of 
p of radius r, written B r (p). When G is not a triangle group, if r < r^ 
then for g G G, g(B r (p)) fl B r (p) ^ if and only if g is in the stabilizer 
of p , jlj. Consequently, the projection of B r (p) is exactly fc-to-l in the 
complement of p. For larger values of r, the projection of B r (p) can be 
considerably more complicated. 

Let 7 : (— oo, oo) — » S be a geodesic parameterized by arc length. 
A lift 7 of 7 to EI has endpoints 7+ and 7_ in the extended real line 
M = MU {00}, representing respectively the limits at infinity and 
minus infinity. If the domain of 7 is restricted to [0, 00) then we shall 
refer to it as a geodesic ray. A lift of the ray 7 has the single endpoint 
7 + at infinity. Henceforth, all geodesies shall be parameterized by arc 
length. 



2.1.1. The distinguished ray X. It is necessary to choose a distinguished, 
simple geodesic ray A on S with initial point A(0) = p, which is used to 
catalogue the excursions of a geodesic on S into a cone neighborhood 
of p . The two cases, S non-compact and S compact, are considered 
separately. 
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If S is non-compact then it has a cusp end, also called a puncture, 
and there is at least one simple geodesic rays with initial point p that 
eventually go out the cusp end. Choose such a ray and call it A. 

Now suppose S is compact. We have assumed that G is not a triangle 
group and therefore, on S there is a simple closed geodesic (5 that does 
not contain p . We shall choose A to be a simple geodesic ray that 
twirls asymptotically into f3. It can be constructed as follows. Let 5 
be a minimal length, and therefore simple, geodesic arc with its initial 
point at p and its terminal point on f3. Let (3 : M — » EI be a full 
connected preimage of (5. Then there are lifts p of p and 5 of 5, so that 
S has its initial point at p and its terminal point on (3. Let A be the 
geodesic ray from p to j3 + , one of the endpoints of f3. Then one easily 
shows that for g G G, either g(X) fl A = or else the intersection is the 
point set {p} and g is in the stabilizer of p . The ray A thus projects 
to a simple geodesic ray on S with initial point p . 

Observe that, given a sequence of values tj — > oo, in the non-compact 
case X(tj) has no limit points in S, while in the compact case the 
sequence \(tj) will limit exactly at points on the geodesic (5. 

2.1.2. Normalization of the covering group. For the sake of clarity and 
computational simplicity we shall choose a nicely normalized Fuchsian 
group G defining the orbifold S . First, let the complex number i G H 
be our point p , covering the cone point p . Then the transformation 



T k (0 = 



cos f C + sin f 



- sin | C + cos f 

generates the stabilizer of % in G. 

We may further normalize G so that the geodesic A has a lift Ao 
with initial point i and terminal point G M. Then Ai = Tk(Xo) 
and A_i = T^T 1 (Ao) are two other lifts of A beginning at i. The point 
tan | = a k is the endpoint at infinity of Ai and — a k is the endpoint of 
A-i- 

We set some notation. Define the sets / = [(— a k , 0) x (0, oo)] U 
[(0, a k ) x (-oo, 0)] and J = [{-a k , 0) x (a k , oo)] U [(0, a k ) x (-oo, -a k )\. 
We shall use the shorthand B r for the disc B r (i) and B r for the disc 
B r (p). 

2.2. Excursions. 

2.2.1. The definition. Given a cone neighborhood B r of p and a geo- 
desic ray 7 on S , we are concerned with the "excursions" of 7 into 
B r . We saw in the introduction that these can be regarded as certain 
arcs of 7, although their definition is complicated by the self-overlap of 
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B r for larger values of r. Hence, we would like to distinguish the fine 
structure of the excursions as 7 intersects different overlapping pieces 
of B r . 

The above considerations dictate the precise, but non- intrinsic, def- 
inition of an excursion which follows. Suppose r > 0. Define an r- 
excursion e of a geodesic 7 into the cone neighborhood B(r) to be a 
lift 7 of 7 to H with (7+, 7_) G / so that 7 n B r ^ 0. 

The r-excursion e = 7 is called an approximating r-excursion if 
(7 + ,7_) G J. More generally, call e an approximating excursion if e is 
an approximating r-excursion for some r > 0. In Section 0221 we shall 
prove that every approximating excursion is an r-excursion for some 
r < r& = sinh - (cot j). 

Let 5(7, r) be the set of r-excursions of 7 and let ^(7) be the set of 
approximating excursions of 7. 

2.2.2. Ordering the excursions along 7 . To each r-excursion e = 7 
there is a single intersection of the geodesies 7 with the ray Ao- We may 
therefore associate to e the unique real parameter t e , for which j(t e ) G 
Ao- Given a positive real value r and a geodesic 7 the identification 
e — > t e defines a function ^ : E(j,r) — > (—00,00). Let E + {j,r) be 
the subset of e G E{^,r) so that ^(e) > 0. ^ a ^ so induces a map on 
A(7) and we similarly define A + ( / y). The range of if) shall be called the 
excursion parameters, written {t e }. We then have 

Proposition 1. If e and e' are two distinct r-excursions along 7 , then 
^(e)^^(e'). 

Proof. Write e = 7, e = 7' and suppose t e = t e /. Since 7 is distinct 
from 7', there is a non-trivial g £ G so that ^(7) = 7' and ^(7(t e )) = 
7'(t e /). It follows that (?(Ao) fl Ao 7^ 0. This is only possible if g — Tj? 
for some n = 1, . . . , k — 1. But then because we have (7 + ,7_) G I, 
applying the transformation gives ((7(7+), ^(7-)) = (7+,7^_) ^ which 
is impossible. 

□ 

Define an order on E(j, r) by stipulating that e < e' if ^(e) < ^(e'). 
This ordering does not depend on the parametrization of 7. The next 
proposition shows that this ordering of E is very well behaved. 

Proposition 2. Given a geodesic 7 , which is distinct from the simple 
closed geodesic (3 used to define X, suppose that the set of r-excursions 
E + = E + ('-f,r) is infinite. Then there is a unique map from IN onto 
E + , making E + into a sequence {ej}'jL 1 so that Cj < et if and only if 
j < k. Furthermore, lim ip(ej) = 00. 

j— >oo 
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Proof. As G is discrete, E + must be countable. Therefore, it will suffice 
to show that for any sequence {e J }°^ 1 in E + , which is ordered as in the 
proposition, lim^oo ip(ej) = oo. We argue by contradiction. Suppose 
there is such a sequence of r-excursions for which the excursion param- 
eters, which we shall write ip{ej) = tj, do not diverge to oo. By passing 
to a subsequence, of the same name, we may suppose tj — > t* < oo. 

Write Cj = jj. Choose a lift 7 of 7. For each j G IN there is a transfor- 
mation gj G G with gji'jj) = 7. These transformations must be distinct, 
since they respect the parameterization; in particular, <7j(7j(0)) = 7(0). 
Then gj(B r ) fl 7 7^ and j(tj) = gj(Xo) fl 7. We take the liberty to de- 
fine, hopefully without confusion, gj(Xo) = Xj- The rays Xj intersecting 
7 in a sequence of points j(tj) converging to j(t*). 

Xj has endpoints gj(i) G EI and bj G IR. Since the interiors of the 
rays are all disjoint and the sequence {tj} is increasing, the sequence 
{bj} can limit at, at most two points in G IR, on opposite sides of 7 . 
By restricting to one side, we may further stipulate that the sequence 
{ej} was chosen so that bj — > b* G IR. 

The initial point gj(i) of the ray Xj is the center of the translate 
gj{B r ) of the disc B r . Since G is discrete, the points gj(i) will limit 
on IR and not in H. Consequently the discs gj{B r ) have Euclidean 
radii going to zero and they also will limit on IR. Recall that for each 
j, gj(B r ) fl 7 7^ 0. Therefore the points gj(i) all lie in the radius r 
neighborhood of 7 and must limit at an endpoint of 7, denoted by %. 

If the rays Xj accumulate in H, it must be at a lift of /3. The limit 
point 6* cannot be an endpoint of 7 distinct from 7*, for then the Aj 
would accumulate on 7 , contradicting the hypothesis. Thus if b* is an 
endpoint of 7 then 6* =7*. 

Putting this all together: we have the sequence of geodesic rays Xj, 
with the endpoints on IR converging to b* and the endpoints in EI con- 
verging to 7*. In order for this to happen the sequence of intersections 
Xj fl 7 = j(tj) must converge to 7*, contrary to the assumption that 
the excursion parameters ip(ej) = tj, do not diverge to 00. 

□ 

If the ray 7 contains infinitely many approximating excursions then 
they are ordered as a subsequence of the r fc -excursions. We write 
{e*}°^ x for the sequence of approximating excursions in ^(7). 

2.2.3. The depth of an excursion. Given a geodesic a G H, define d(a) 
to be the distance from a to the point i. Define the depth of the r- 
excursion e = 7 of 7 , to be the value d(e) = ^(7). For a pair of points 
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(x,y) G /, let a be the geodesic with endpoints (a+, a_) = (x, y) and 
define D(x,y) = d(a). 

3. The application of Ergodic Theory 

3.1. The geodesic flow. The unit tangent bundle TiH over EI may 
be identified with HxS 1 . Define the measure fl = A x 9, where A and 
9 are respectively, area measure on H and lebesgue measure on S l . The 
measure fl is invariant under the geodesic flow & on TiH. 

A full measure subset of TiH can be modeled by triples (ip, C, t) G IR 3 
where if) ^ (. Let a be a geodesic in H parameterized so that a(0) is the 
Euclidean midpoint of the semi-circle a(M). If if> = a + and £ = a_, 
then (if),(,t) corresponds to the unit tangent vector a(t) G TiH. It is 
clear that the correspondence defines an injection onto the subset of 
non- vertical vectors in TiH. In these coordinates the invariant measure 

u has the form u = — —dipd(dt, up to a constant multiple. 

The measure fl projects to a measure /i on TiS, which is invariant 
under the the geodesic flow G* on TiS. It is well known that when 
S is a finite area surface, & acts ergodically with respect to \i. Every 
orbifold has a finite branched cover that is a hyperbolic surface. Thus 
ergodicity of the flow on the surface implies the ergodicity of the flow 
on the orbifold. See [T], [TU] and [TT] for details. 

3.2. Convergence of the distribution. Every v q G TiS uniquely de- 
termines a geodesic (or a geodesic ray) 7 where 7(0) = q and 7(0) = v q , 
the unit vector tangent to 7 at q. Conversely, every geodesic determines 
a unit tangent vector at the parameter t = 0. We say that a property 
holds for almost all geodesies on S if there is a set A C TxS, of full 
/x-measure, so that the property holds for every geodesic determined 
by a vector in A, as above. 

There are two classes of subsets of / and J that play particularly 
important roles in what follows. For < z < 00 set Q(z) = {(x,y) G 
i" I D(x, y) < z} and for < z < r k set fl*(z) = {(x, y) G J\ D(x, y) < 
z}. As a consequence of Corollary |2l £l*(z) = J for z > r^. At times we 
shall use the letters ip and ( in place of x and y. Then for appropriate 
values of z define the integrals 
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After the next theorem, the proofs of Theorems Q and El become 
computations in hyperbolic geometry and plane integration. 

Theorem 3. For almost all geodesic rays 7 on S , for r > and 

< z < r, distk{r, z)^) converges to distk{r, z), where 

dist k (r,z) = ^j. 

For almost all geodesic rays 7 on S and for < z < r&, dist* k (z)(^) 
converges to distl(z), where 

The remainder of Section El shall be devoted to the proof of this 
theorem. The idea is to produce a measure transverse to the geodesic 
flow, on which there is an ergodic action related to the sequence of 
excursions along a generic geodesic on S . The distributions can be 
expressed as limiting sums with respect to this action which then, using 
the Birkhoff Ergodic Theorem, become the above integrals. 

3.3. A cross-section of G l over T% A. Let T\ A denote the subset of the 
unit tangent bundle over the geodesic A, which consists of the vectors 
v q with q G A \ {p}. We shall define subsets of 7\A associated with 
r-excursions and separately , with approximating excursions. Most of 
the work shall take place in H, on the preimage Ti\ of 7\A. 

Given r > define the set of vectors V C T1A0, where v q G V if there 
is a geodesic a in EI with 

(1) (ar + ,a_) G /, 

(2) a(0) = q G A \ {i} and d(0) = v q , 

(3) a(lB r ^ 0, and 

(4) the geodesic ray noa on S contains infinitely many r-excursions. 

Similarly, define /* C TiA , where v q G /* if there is a geodesic a in 
H with 

(1) (a+,a_) G J, 

(2) a(0) = q G A \ {«} and d(0) = v q , and 

(3) the geodesic ray 7r o a on S contains infinitely many approxi- 
mating excursions. 

In order to treat the various cases simultaneously, we adopt the con- 
vention that r denotes either the parameter r > or the symbol *. 
The set of vectors l T projects to a set of vectors l T C 7\ A. Essentially, 
l T consists of the tangent vectors of the form j(t e ), associated to excur- 
sions e, with the caveat that the ray 7 contains infinitely many such 
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excursions. Since lis simple, the projection 7T* : Ti\ — > TiA is bijective 
and consequently we have the following. 

Proposition 3. The projection 7r* : V — > Z T is a bijection. 

We let *-excursion be another term for approximating excursion. 
The next proposition asserts that l T is a cross-section for the geodesic 
flow on S . 

Proposition 4. For each value r, there is a full measure, G t -invariant 
subset M T C TiS so that V is a cross-section for the flow G l acting 
on M T . In particular, for each v G M T , there is a geodesic ray 7, 
possessing an infinite sequence of r-excursion parameters {t e }, so that 
7(0) = v and G l (v) G V fort > if and only ift = t e . for some j G IN. 

Proof. We argue when r — r. Let M r be the set of vectors v G T\S 
so that a geodesic ray 7 with 7(0) = v contains infinitely many r- 
excursions, none of which intersect the cone point p. It follows from 
Proposition El that M r is G*-invariant. In particular, l r C M r . 

In order to see that M r has full /^-measure, choose a value u so that 
< u < rfc and u < r. As a consequence of the Poincare Recurrence 
Theorem [7j, there is a set of full measure N u C XiS 1 , with the property 
that for each v G N u the geodesic 7 with 7(0) = v returns infinitely 
often to B u . It is proved in Corollary El of Section 14.2.21 that every 
intersection of a geodesic 7 with B u produces a -u-excursion. Therefore, 
such a geodesic 7 must contain an infinite sequence of w-excursions. 
The subset iV* consisting of v G N u so that the G*-orbit of v does not 
pass thru p, has full measure in N u . As u < r, N* C M r , showing that 
M r has full measure in T\S. 

Given v G M r , let 7 be the geodesic with 7(0) = v. Then 7 contains 
an infinite sequence of r-excursions {ej} with excursion parameters 
ip{ej) = t e .. As before write ej = jj with ((7^)+, (jj)-) G I. Thus for 
each j G W, (v) = ^{t e .) G / r . 

If for some t > and t> G Af r , G l (v) G / r , then there is a lift v of 
v and 7 of 7 so that j(t) — v G I\Aq, (7 + ,7_) G / and 7 fl 5 r 7^ 0. 
Therefore e = 7 is an r-excursion. The argument for approximating 
excursions is similar. 

□ 

3.4. The first return map. Let F T : V — ► V be the first return map 
under the geodesic flow G l : M T — > M r . In other words, given o £ l T let 
7 be a geodesic in 5* with 7(0) = v. Then F T (v) = j(t), where t > is 
the first value with ^(t) G V . The invariant measure \i for the geodesic 
flow on S induces an F r -invariant measure v T on l T for which F T is 
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ergodic, (UI21- Under the identification of V with l T , F T lifts to a map 
F T : l T —>■ l T which is invariant and ergodic with respect to the lifted 
measure u T . Furthermore, this construction produces an isomorphism 
between the dynamical systems (Z T , F T , v T ) and (7 T , F T , u T ). 

It will be easier to work with F T in a form associated to the (ijj, £, t) 
coordinates. Let Q T denote Q(r) if r = r and J is r = *. There is a 
map B T : l T —> fi T , taking v q to the endpoints (a+, a_) of the geodesic 
a with a(0) = q and d(0) = v q . B T is injective in the complement 
of the subset of l T of periodic points of the map F T . Since the set of 
periodic points is countable, B T is injective on a set of full //-measure. 
Henceforth, we shall remove the set of periodic points from V and V ' . 

If we represent v q = (ip q , ( q , t q ) in coordinates, then B T becomes 
the projection B T (if} q ,Q,t q ) = (ip q ,Cq)i onto the first two coordinates. 
Let u) T = B T (l T ). and define the transformations R T : uj t — > uj t by 
R T o B T = B T o F T . Note that uj t is precisely the subset of Q T so that if 
a is a geodesic with (a + , a_) = (V 1 , C) £ the geodesic ray a projects 
to a geodesic ray on S that contains infinitely many r-excursions or 
approximating excursions, not thru p, as appropriate. 

Arguing as in jTT] , we can use a theorem of Ambrose 3] to show that 
the geodesic flow can be built up, in a nice way, from the first-return 
map on the cross section. This theorem allows us to conclude that 
the G*-invariant measure /x on S has the form /j = v T x dt in local 
coordinates, where dt represents lebesgue measure on the flow lines. 
Then // = v T X dt pulls back to the G*-invariant measure jl = v T x dt 
on T]H. It follows that in (/0, ^) coordinates, v T = T^z^p li^^C is the 
induced invariant measure for i? T . We shall write v T = v since, in these 
coordinates, the measure is independent of r. 

Since M T is full measure in T\S and G*-invariant, the cross-section 
l T is of full measure in 7\A. It follows, in particular, that c<j r is full 
measure in Q T . In light of this, the transformation R T : Q T —>■ Q T is 
well defined up to sets of measure zero, is i? T -conjugate to F T , and has 
invariant measure u T . 

It is also important to note that, as a consequence of Proposition HI 
for almost all v G TiS, there are infinitely many values tj — > oo so that 
the geodesic 7 with 7(0) = v has j(tj) G / r . Therefore, for almost 
all v G TxS the geodesic 7 with 7(0) = v has a lift with its endpoints 
in u T . In other words, almost every geodesic on S has a lift with its 
endpoints in u T . 

3.5. The proof of Theorem^ As usual, we shall prove the theorem 
for r-excursions. The details for approximating excursions are similar. 
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First observe that for (i[) , Co) G Q(r) we have the equality 

(2) lim \l<j<n,Do R^ , Co) < z} = 

n^oo Tl 

1 H 

(3) lim -r( X[M o5)(i?(^C )), 

3=1 

The measure i>* = jhzV is an i?-invariant probability measure on Q(r), 

with respect to which R is ergodic. By Birkhoff 's Ergodic Theorem and 
the fact that u r has full z>*-measure in f2(r), for almost all (V'ojCo) £ 
fl(r) the limit is equal to 

1 f A(Y) 

Almost all geodesies 7 on S* have lifts with endpoints in u T . Let 7 be 
one of these geodesies. Then there is an infinite sequence tj = jj, j = 
0, 1, ... , of r-excursions along the ray 7 , with associated excursion 
parameters tj = 4>(ej). Let Vj = j(tj) G V . By Proposition HI F T (vj) = 
Vj + i and therefore F^(v ) = Vj. The lift of Vj to A is Vj = 7j(£j) and, 
as above, F^(vo) = Vj. 

Write 7 = 7 . To each j G IN there is a unique gj G G so that 

7j = %(7)- Th en in (V>,C>*) coordinates, = ((-%•)+, (7j)-> s i) = 
(<7j(7 + ), ^(7_), Sj), for some values s 3 - that are not related to the tj. 
Then if we set (?/>o, Co) = (7+, 7-) G c<j t we get 

(4) i#(Vb, Co) = 5r o o b;\%, T _) = fi r O i^ ) = 

5 T (u,) = J B r (^(^ ),^(Co),s j ) = (gj(ip ),9j(Co))- 
Let Xy denote the characteristic function of the set Y. Turning to 
the distribution we have 

dist k (r, z)(j) = lim 1 1 < j < n, dfe) < z} = 

n^oo n 

lim \l<j<n, dtfj) < z} 

lim I 1 < j < n,D( 9j ^ + ), gj (^)) <z} = 

n^oo Tl 

(5) lim I 1 < j < n, D(g^ ) } 9j (( )) < z}. 

n— »oc n 

By the sequence of equalities (HJ), the limit (J3j) is equal to the limit 
(J2J). Therefore, we may conclude that for almost all geodesies 7 on S , 
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dist(r, z)(j) converges to dist(r, z) = ^14, which proves Theorem 01 for 
r-excursions. □ 



4. The computation of A and A* 

4.1. The values of A and A*. In light of Theorem the proofs 
of Theorems ^ and 121 will be completed by computing the values of 
A and A*. Recall that Vk = sinh _1 (cot |), 5k = sinh _1 (cot ^) and 

/ n 1— xtan-f 



Theorem 4. For z > we have 



A(z) 



A*(z) 



^ sinh z if z < r\ 

sinh z tan -1 ( ^ ) + log(sin | cosh z) if z > r^, 

sinh z if z < 5k 

sinh z tan -1 (ipk (sinh z)) + log (2 cosh z sin | cos |) if cifc < z < 



log(2cosf) if « > r fc . 

Mien = 3 or 4, 5k < and only the second two cases occur for the 
value of A*(z). 

The rest of this section is devoted to the proof of this theorem. 
4.2. Hyperbolic geometric considerations. 

4.2.1. The geodesic tangent to a disc. In order to turn A and A* into 
easily computable double integrals, we shall derive several formulae 
that will be of use in determining the limits of integration. Given 
distinct points a, b G IR, let ah denote the geodesic a in HI with a + = a 
and o_ = b. 

Theorem 5. Given x > 0, the point w G [—-,x) for which the geodesic 
xw, is tangent to the hyperbolic disc B p is given by the formula 

x sinh p - 1 

w = Wpix) = — . 

x + sinh p 

Similarly, given x < the point w G (x, —-] for which the geodesic 
xw, is tangent to the hyperbolic disc B p is given by the formula w = 
-W p (-x). 
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Let A denote the unit disc model for the hyperbolic plane, which 
we envision as sitting in the complex plane. Write <9A for the unit 
circle, which is its boundary at infinity. The proof of Theorem [5] in- 
volves computations in A and the transformation of those results to H. 
The following lemma can be proved by an easy computation using the 
hyperbolic metric in A. 

Lemma 1. In A, the circle of hyperbolic radius p and center has 
Euclidean radius r = tanh | 

If z, £ G dA then, as in the previous case, we let z£ be the hyperbolic 
geodesic in A with endpoints z, £. Similarly, for z G dA let zO be the 
geodesic with endpoint z passing thru the origin. This geodesic is a 
Euclidean straight line. We use the same notation B p (0) to denote the 
hyperbolic disc of center and radius p in A. To avoid confusion we 
shall continue to write B p for B p (i). 

Lemma 2. Suppose that z,£ G dA are two points so that the geodesic 
zl; is tangent to the the hyperbolic disc B p (0). Then \z — £| = Isechp. 

Proof. There is no loss of generality in supposing that z = x + iy with 
x,y > and £ = z. Then we have \z — £| = 2y and it remains to 
show that y = sechp. Let C be the circle perpendicular to dA whose 
intersection with A is the geodesic z£. Since C is invariant under 
the reflection z — > = which fixes dA, it crosses the real axis in points 
< r < 1 and -. By Lemma EJ since C is tangent to the circle of radius 
p, r = tanh |. 

As a consequence of z lying on C, it satisfies the equation \z — + 
r)| 2 = — r)] 2 . Also, since z G dA, \z\ = 1. Solving these equations 
simultaneously gives y = Setting r = tanh | and simplifying, 

gives the result. □ 

Lemma 3. Suppose z is a point in the upper half of the unit circle 
dA. Let £ G dA be the point which lies to the right of the line zO so 
that zt; is tangent to the disc B p (0). Then £ = z{c — i\J\ — c 2 ) where 
c = 1 — 2sech 2 p. The above also holds for z = 1, where £ is then a point 
in the lower half of dA. 

Proof. From the previous lemma we have \z — £| = 2sechp, which sim- 
plifies to Rez£ = 1 — 2sech p. For simplicity write c = 1 — 2sech p. Note 
that as p increases on the interval [0, oo), c increases on [—1,1). 

The point £ may be written in the form £ = rjz, for a point r\ of 
modulus one. Let r] = u + iv . One see by a simple computation that 
u = Rezl; = c and then w = ±vl — c 2 . With the choice of the minus 
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sign the points rj = c — iy/1 — c 2 fill out the lower half of the unit 
circle. As a function of p, r\ moves counter-clockwise through its arc as 
p increases from zero to infinity. Thus, treating the point £ = rjz as a 
function of p, we see that it moves counter-clockwise around the circle 
from —z to z as p increases from zero to infinity. It follows that for all 
values of p, the point £ = rjz lies to the right of the line zO, or in the 
lower half-plane if z — 1. With the choice of a plus sign, rj fills out the 
upper half of the circle and the points rjz lie to the left of zO, or in the 
upper half-plane if z = 1. That completes the proof. □ 

Proof of Theorem^ The Mobius transformation g(z) = ~ % Z Z X/ maps A 
iso metrically to H, taking the disc B p (Q) to the disc B p . g also defines 
a one-to-one correspondence between dA and IR, where 1 goes to 0, 
-1 goes to oo and the upper and lower semi-circles respectively, are 
mapped to the positive and negative real axes. 

As in the previous argument let 77 = c — — c 2 , where c = 
1 — 2sech 2 p, and define the transformation W p (x) = g{r]g^ l {x)). Sim- 
plifying we get 



which shows that the formula for W p is the one asserted in the theorem. 
Now we need to see that W p does what it is claimed to do. 

If x > then g~ 1 (x) = z lies in the upper-half of dA or is equal to 1. 
Then £ = rjz is the point of dA for which the geodesic z^ satisfies the 
hypotheses of the previous lemma. The image of this geodesic, g(z£), is 
the hyperbolic geodesic in EI which is tangent to g(B p (0)) = B p and has 
endpoints x and W p (x). Recall from the proof of Lemma El that £ = rjz 
lies in the counter-clockwise arc of dA between —z and z. Since the 
image under g of this arc is the interval [—-,x), the proof is complete 
for x > 0. 

Now suppose x < 0. Consider the isometry h(z) = —z of H. h 
maps the geodesic tangent to B p , with endpoints — x and W p (—x) G 
[-, —x) to the geodesic tangent to B p with endpoints x and — W p (— x) G 
(x, —-]. Thus — W p (— x) is as asserted in the proposition. □ 




A further computation gives 
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4.2.2. Consequences of the tangency computations. 

Corollary 1. Given x, z > 0, D(x,y) < z for y G [—-,x) if and only 

ifye[-i,W z (x)}. 

Proof. For fixed x > and for y G [—~,x), define H(y) = D(x,y). 
Then H(y) is the value p for which the geodesic xy is tangent to the disc 
B p . We have shown in the proofs of Lemma El and Theorem 03 that H 
is an increasing function of y, mapping [—-, x) onto [0, oo). Also recall 

that, for x, y as above, xy is tangent to B z precisely when y = W z (x). 
It follows that D(x, y) = H(y) < z if and only if y G [— ^, W^(a;)]. □ 

Corollary 2. When k > 5, every z-excursion with z < 5 k is an ap- 
proximating excursion. Moreover, every approximating excursion is a 
z-excursion for some z <r a . 

Proof. First observe that the equation W z {a k ) = —a k has the solution 

z = smh" 1 ^ 1 ^)) = sinh"^ 1 "^* ) = sinh _1 (cot ^) = 5 k . 

This is positive if and only if k > 5, which we henceforth suppose to 
be the case. For x, z > 0, W z (x) is an increasing function in both 
variables. Therefore for z < 5k and x G (0,a&), W z (x) < W z {(ik) < 
W Sk {a k ) = -a k . 

Let e = 7 be a z-excursion for some z < 5 k . Then x = 7 + G (0, a k ) 
and, as a consequence of the above, 7_ = W^(x) < — a k . This says 
that e is an approximating excursion, proving the first assertion of the 
corollary. 

Now observe that z = D(0, —a k ) is the value for which W 2 (0) = — a k . 
Solving for z we get z = sinh(cot(|)) = r k . Thus the geodesic a with 
endpoints and a k is tangent to the disc B rk . 

Let e = 7 be an approximating excursion. Then 7+ G (0, a k ) and 
7_ G [— -, — a k ). It follows that 7 lies above a in EI and therefore, 
7 fl B Tk 7^ 0. The second part of the corollary is proved. 

□ 

Corollary 3. Suppose 7 is a geodesic on S and 7 D B rk = e 7^ 0. Then 
there is a lift 70/7 and a lift e of e contained in^ , so that e\~\B r . k 7^ 
and e = 7 is a z-excursion for some z < r k . 

Proof. Choose a lift 7 so that £nB rk 7^ and 7+ G (0, a k ). Without loss 
of generality we suppose that 7_ < 7+. Since W rk (a k ) = 0, the geodesic 
a with endpoints and a k is tangent to the disc B n . If 7_ > then 
7 lies entirely below a and 7 fl B Tk 7^ 0, which is impossible. Therefore 
7_ < and e = 7 is a z-excursion for some z < r k . □ 



The distribution of geodesic excursions 



17 



4.3. The proof of Theorem 01 We shall write the integrals as sums 
of iterated double integrals. Once in this form, their actual evaluation 
is a elementary. 

4.3.1. The computation of A(z). Consider fixed values z > and x G 
(0,a h ). If W z (x) < then, by Corollary □ D(x,y) < z for y G [-~,0) 
precisely when y G [— -, W^x)], whereas if W z (x) > then D(x, y) < z 
for all y G [— -, 0). As a result of this observation, for fixed z > 0, there 
are two cases to be considered. 

The first case is when z < r k . This condition is equivalent to a k < 
-r^r - or W z (a k ) < 0. Since W z (x) is an increasing function of x, this 
occurs if and only if W z (x) < for all x G (0, a^). Then, as observed 
above, for x G (0,a fc ), D(x,y) < z for y G [— -,0) if and only if y G 
[ — , W z (x)]. By symmetry, we have that for x G (—0,^,0), D(x,y) < z 
for y G (0, — -] if and only if y G [— W z {— x), — -]. Consequently, when 
z < ffc we get 

/•a* /.W»(aO j /•- a fc /— i 1 

A(z) = / r^dydx + / / -^dydx 

Jo J-\ \x ~ VY Jo J-w z (-x) [x - y) 2 

In the second case z > r k and we will have W z (x) = for x = 
(0, afc). Then the set of y values for which D(x,y) < z will 



sinh z 

take two different forms. For x G (0, in>i ), we have D(x,y) < z for 
y G [-J,0) if and only if y G W z (a;)]; whereas, for x G [^rj,0) 
we get D(x,y) < z for all y G [— ~,0). Similar results again hold, by 
symmetry, for x G {—a, 0). Then for 2 > r k we get 



o 







1 



dydx + / —dydx 



4.3.2. T/ie computation of A*(z). The same considerations and cases 
come into play in the computation of A* {£). In addition there is a 
third special case, which is a consequence of the fact that every ap- 
proximating excursion is an r-excursion for some r < r k . The three 
cases depend on the position of the points W z {x) with respect to — a k , 
rather than with respect to 0, as in the earlier computations. The cases 
are also described by the position of the particular point x z , for which 
W z (x z ) = -a k . 

In the first case z < 5 k . As observed in the proof of Corollary |2J 
this is equivalent to W z (x) < —a k for all x G (0,a k ). Since 5 k < 
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for k = 3,4, this case does not arise for those values of k. Then, 
by Corollary [H for all x G (0,a k ), D(x,y) < z for y G [— — a k ) if 
and only if y G [—-, W z (x)], and by symmetry, we have that for all x G 
(-a k , 0), D(x, y) < z for y G (a fe , -±] if and only if y G [-W*(-ar), --]. 
This is the same situation we had computing the integral in the earlier 
first case. Therefore for z < 6 k , 

A(z) = A*(z) = sinh z. 

k 

In the second case 5 k < z < r k , where we let S k = for k = 3,4. 
Consider the equation W z (x) = —a k , whose solution is a point x z = 
(p k (smh z) . x z is a decreasing function of z. We've seen that xs k = a k 
and x rk = 0. Therefore, for 5 k < z < r k we have W z (x) = —a k 
for the point x = x z G (0, a k ). As in case two in the computation of 
A, the set of y values for which D(x,y) < z will take two different 
forms. For x G (0,x z ), we have D(x,y) < z for y G [— -, — a k ) if 
and only if y G [— -,W z (:c)] and for x G [x z ,0), we get D(x,y) < z 
for all ?/ ^ [ — \i —^k)- Similar results again hold, by symmetry, for 
x G (— a k , 0). It follows that for S k < z < r k 

A " (2) = X L i^w dydx + L L j^w dydx + 



o 



i , , r~ Xz r~ i 



dydx + / —dydx 



-w z (~x) (x - y) 2 J-a k Jo (x - y) 1 

The last case to consider is when z > r k . By Corollary [21 £l*{z) = J 
for all z > r k . Therefore, for all x G (0,a k ), D(x,y) < z for all y G 
[ — , —a k ). Consequently, in this last case we have, for z > r k 

ra h r-a k ^ /■- o* /■- 1 ^ 

A*(z) = / / r^dydx + / r^dydx 



x-y) 2 ' Jo Ja k (x-y) 

5. Depth as a function of area 

We have been dealing with excursion depth as expressed in terms of 
the radius of the largest cone neighborhood disjoint from the excursion. 
In jS], where we looked at excursion into a cusp, it was just not possible 
to express the depth of an excursion as the distance to the cusp end, 
since the cusp is off at infinity. Instead we considered the area of the 
largest neighborhood of the cusp disjoint from the excursion. Suppose 
we were to take the same approach here. 

In the hyperbolic plane EI the area of a hyperbolic disc of radius r is 
27r(coshr — 1) [3]. Let G be a Fuchsian group representing S and let 
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Ti : EI — > M/G = S denote the natural projection map. Suppose p is 
the cone point on S and p is a preimage in HI. The actual area of a disc 
-B r (p) is at least 1/fc times the area of its preimage B r (p), with equality 
if the projection is precisely k-to-1. In the case of equality, if A is the 
area of B r (p), then B r (p) has area kA = 27r(coshr — 1). Even if there 
is overlap and the projection is not k-to-1, in order to be consistent 
with |5], we define the A-depth of an excursion e, written A(e), to be 
^ L (cosh((i(e)) — 1). Define an /^-excursion e to be a usual r-excursion 
where r = cosh _1 (^:i? +1), the radius of the disc of area R about p . 
The definition of an approximating excursion is unchanged. 

Suppose {e, } is a sequence of i?-excursions along a geodesic 7 on S . 
Then define 

Adist k (R, Z)(j) = lim \l<j<n, A{e j ) < Z}. 

As a consequence of the definition of A(e), this is equal to 

1 k 
lim 1 1 < j < n, d(e,j) < cosh _1 (— Z + 1)}. 

?woo n 2ir 

For almost all geodesies 7 on S this last limit was shown to converge 
to 

+ D.cosh-^Az + 1)) = ^||f^ . 

For R < ^(cosh(rfc) — 1) and Z < R, the above takes the form 



2tt\2 _ ( 2n\2 
k I \ k ) 



(R + 



2tt \2 _ I 2vr\2 ' 
k > \ k > 



which is almost linear in Z and which approaches Z/R as k — * 00. 
This last expression in the variable Z is the distribution function for 
an i?-excursions along a geodesic with respect to a cusped end of a 
surface, when R < 2, [Hj. Note also that as k — > 00, Rk approaches 2. 

For approximating excursion, the corresponding distribution in terms 
of area has the form 



21og(2cosf) 



for Z < ^(cosh(5fc) — 1). Again as k — > 00 the distributions limit 
at Z/21og2, the value of the distribution of approximating excursions 
alone geodesies out a cusped end when Z < 1. 
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Clearly, the other cases look rather nasty. Nevertheless, one can 
verify that they limit at the corresponding distributions on surfaces 
with cusped ends. 
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